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ABSTRACT 

We discuss massive outflows in galaxy bulges, particularly ones driven by accretion 
episodes where the central supcrmassive black hole reaches the Eddington limit. We 
show that the quasar radiation field Compton-cools the wind shock until this reaches 
distances ^ 1 kpc from the black hole, but becomes too dilute to do this at larger radii. 
Radiative processes cannot cool the shocked gas within the flow time at any radius. 
Outflows are therefore momentum-driven at small radii (as required to explain the 
M — a relation). At large radii they are energy-driven, contrary to recent claims. 

We solve analytically the motion of an energy-driven shell after the central source 
has turned off. This shows that the thermal energy in the shocked wind can drive 
further expansion for a time ~ 10 times longer than the active time of the central 
source. Outflows observed at large radii with no active central source probably result 
from an earlier short (few Myr) active phase of this source. 

Key Vifords: galaxies:evolution - quasars: general - black hole physics - accretion 
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1 INTRODUCTION 

The large-scale structure of galaxies often has surprisingly 
close connections to properties of their nuclei. The M — a 
relation between supermassive black hole (SMBH) mass M 
and bulge velocity dispersion a is the most striking of these. 
Similar relations hold between black hole and galaxy stellar 
bulge mass Ait,, and between the mass of nuclear star clusters 
and a in galaxies where there is no strong evidence for the 
presence of an SMBH (ct < 150 kms~^). 

Massive gas outflows driven by the central object of- 
fer a way of connecting these apparently disparate scales. A 
fast wind from the nucleus collides with the host galaxy's 
interstellar medium, driving a reverse shock into the wind, 
and a forward shock into the ISM. This shock pattern moves 
outwards at a speed mainly determined by whether or not 
the reverse shock cools on a time short compared with the 
outflow timescale (R/R) or not (cf Dyson & Williams, 1997; 
Lamers & Cassinelli, 1997). In the first case (efficient cool- 
ing) , only the ram pressure of the original outflow is commu- 
nicated to the ambient medium. This is a momentum-dnven 
flow. In the second case (inefficient cooling) the full energy 
of the fast wind is communicated to the ambient medium 
through its thermal expansion after the shock. This is an 
energy-driven flow, which expands at higher speed and so 
can have a much larger effect on the bulge of the host galaxy. 

Both types of outflow are important in galaxy forma- 
tion. This paper is mainly concerned with the large-scale 
effects of energy-driven flows. The existence of flows of this 
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type has recently been questioned, so we flrst set the prob- 
lem in context. We show that energy-driven outflows do oc- 
cur, and are ubiquitous on large scales. Solving the outflow 
equations analytically, we give a simple relation between the 
time that the outflow is driven by the central source, and 
the time over which it can be observed as coasting after this 
source turns off. This relation means that observed outflows 
can be used to constrain the past activity of a source. In this 
paper we deal with the case where this source is a quasar, 
which we model as an Eddington-accreting SMBH. Similar 
considerations apply in cases where the driving source is a 
nuclear star cluster. 



2 MOMENTUM OR ENERGY DRIVING? 

The flrst proposal that outflows might relate SMBH and 
galaxy properties was by Silk & Rees (1998, hereafter SR98), 
who considered the effect of an Eddington wind from the 
black hole colliding with the host ISM. Requiring the shock 
pattern to move with the escape velocity, and so presumably 
cutting off accretion to the black hole, they found AI (x a^, 
with an undetermined coefficient of proportionality. Later, 
King (2003; 2005) pointed out that SR98 implicitly assumed 
an energy-driven outflow, whereas Compton cooling in the 
radiation fleld of the active nucleus was likely to produce a 
momentum-driven flow. The condition that this flow should 
be able to escape the immediate vicinity of the black hole, 
and so cut off accretion, predicts a black hole mass 



2 X lO*M0cr|oo 



(1) 
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in good agreement with the observed relation (Ferrarese & 
Merritt, 2000; Gebhardt et al., 2000) (which is itself proba- 
bly an upper limit to the SMBH mass, cf Batcheldor, 2010; 
King, 2010b). Here fg ~ 0.16 is the gas fraction, k the elec- 
tron scattering opacity, and (T200 the velocity dispersion in 
units of 200 kms~^. By contrast, an energy-driven outflow 
as in SR98 would have produced a mass smaller than ([l]) by 
a factor ~ cr/c ~ 10""^ (e.g. King, 2010a). A later application 
of similar ideas (McLaughlin et al., 2006, Nayakshin et al, 
2009) to outflows driven by nuclear star clusters shows that 
these produce an offset M — a relation between total cluster 
mass and velocity dispersion, the offset resulting from the 
fact that star clusters produce roughly 20 times less outflow 
momentum per unit mass compared with an accreting black 
hole. 

There have also been attempts to explain the relation 
between black hole and bulge stellar mass in terms of Ed- 
dington outflows from accreting SMBH. The observed rela- 
tion M ~ W'^Mt, (cf Haring & Rix, 2004) means that this 
an inherently more complex problem than M ~ a, since Mb 
is apparently the small part that remains after some pro- 
cess has almost swept the bulge clear of its orginal baryon 
content. Two recent papers discuss this problem. 

Power et al. (2011) suggest that star formation in a 
galaxy bulge is self-limiting, and this limit largely deter- 
mines the bulge stellar mass Mt- They further suggest that 
an energy-driven outflow from the central black hole clears 
away the remaining gas. This process cannot be totally ef- 
fective: King (2010b) shows that energy-driven outflows are 
Rayleigh-Taylor unstable since the rapid expansion of the 
shocked wind leads to a large density contrast with the am- 
bient medium. Thus a fraction of the gas can still remain 
even after the outflow passes. 

In contrast. Silk & Nusser (2010) assert that energy- 
driven outflows do not occur at all in galaxy bulges. It is easy 
to show that momentum-driven outflows cannot clear the 
remaining gas from the bulge (Silk & Nusser, 2010; Power 
et al., 2011, Appendix). Accordingly Silk & Nusser (2010) 
suggest that star formation must be able to remove it. 



3 SHOCK COOLING 

To decide whether energy-driven outflows exist or not we 
consider an Eddington wind (Mout — M-^dd) from a super- 
massive black hole propagating in an approximately isother- 
mal galaxy bulge, with gas density 



P = 



27rGr2 ' 

The gas mass inside radius R is 



M{R) = 47r / prdr = 



G 



(2) 



(3) 



As we have seen, the important question for the gas motions 
is whether the reverse shock cools. The preshock wind has a 
velocity v~r]c~ 0.1c (King & Pounds, 2003; King, 2010a), 
which implies a (reverse) shock temperature 



Ts 



3 htuh 2 
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1.6 X 10^° K. 
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This gas is clearly too hot to have any bound electrons, so 
the only losses cooling it are Compton and free-free. The 
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Figure 1. Evolution of an energy-driven shock pattern for the 
case (T = 200 kms~^, fg = 10"^ computed numerically from the 
full equation |(TTJ. Top: radius vs time, middle: velocity vs time, 
bottom: velocity vs radius. The curves refer to different initial 
conditions: black solid - Ro = 10 pc, vq = 400 km/s; blue dashed 
- i?0 = 100 pc, VQ = 1000 km/s; red dot-dashed - i?o = 50 
pc, no = 200 km/s. All these solutions converge to the attractor 
JTSj. The vertical dashed line marks the time t = 10^ yr when 
the quasar driving is switched off. All solutions then follow the 
analytic solution I I19I I. 
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Figure 2. Same as Figure 1, but with Ro = 50 pc and vq = 200 
km/s and varying mean gas fractions: blaclc solid - fg = 3 ■ 10~^; 
blue dashed - fg = lO^'^; red dot-dashed - /g = 3 ■ 10~^; green 
dash-triple-dotted - fg = 10~^. 



mass conservation equation for the Eddington outflow gives 
a postshock number density 

N = ^x ~ 1 X 10-='(A/out/A/Edd)M8i?kp'c cm-' 

(5) 

where Ms is the SMBH mass in units of lO^M© and iZkpc is 
the radial distance in kpc. This gives a radiative (free-free) 
cooling time for the shocked gas of 



f,ad =i 2 X lO^Mg-'i^^p, yr. 



(6) 



King (2003, eqn 8) shows that the Compton cooling time of 
this gas in the quasar radiation field is 



- 2 cR^ 
^ 3GM 



V 



(7) 



where 6 ~ 1 is the fractional solid angle of the outflow, 
me,mp are the electron and proton masses, and we have set 
n = 0.1c in the original equation. 

To decide if cooling is effective we compare these 
timescales with the flow timescale for a momentum-driven 
outflow, which is 



R 
R 



5 X 10''i?kpc<T2ooM8 ^'^ yr, 



(8) 



(9) 



(cf. King 2003, eqns 9 & 14). We flnd directly 
= l.9,RvpcCJ2QoM^^''^h. 

t'flow 

We see that Compton cooling is effective only out to 
about R = 1 kpc, (cf Ciotti & Ostriker, 1997) while the 
radiative (free-free) cooling is always far longer than the 
flow time. Silk & Nusser (2010) claim the opposite, but 
appear to have considered the cooling of the ambient gas 
rather than the shocked wind which contains all the en- 
ergy. Their adopted cooling function (Sutherland & Dopita, 
1993) only goes to temperatures 10^ — 10** K, far below the 
shock temperature Ts — 10^" K. We recover the result (King, 
2003; 2005) that in a galaxy bulge an Eddington outflow is 
momentum-driven when very close to the SMBH, but be- 
comes energy-driven outside a typical radius ~ 1 kpc. 

Many galaxies show evidence for massive high-speed 
[v ~ 1000 kms"^) gas outflows on large scales (~ 20 kpc) 
(e.g. Tremonti et al., 2007; Holt et al. 2008). By the argu- 
ments of this Section, these must be energy-driven. Their 
ultimate cause may be starbursts, or AGN activity by the 
central SMBH. However these nuclear phenomena are often 
absent or weak when the outflows are observed. So to un- 
derstand the connection between the observed outflow and 
its original cause we need to know how the outflow coasts 
and ultimately stalls in the absence of driving. 



4 ENERGY-DRIVEN OUTFLOWS 

The equation governing the movement of the shock pattern 
in an energy-driven outflow in an isothermal potential is 
(King, 2005) 



— -t^Edd 



2fga^ 
G 



^r^'r+srrr+^r^ 



+ wh-^R (10) 



with r; ~ 0.1 the accretion efficiency, I/Edd the Eddington 
luminosity of the central black hole, a the velocity disper- 
sion of the ambient medium and fg the gas fraction relative 
to all matter in this medium. The latter quantity may be 
depleted relative to its value fc prevailing when the earlier 
momentum-driven outflow establishes the A'l — a relation 



IT}. Using the expression M — 



m LEdd gives 



J7ca^^ = IIr^R + SRRR+^R'^ 



+ 5aR 



(11) 



This equation has a solution of the form R = v^t with 
J- _ 



2»7C 



(12) 
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[Note that in King (2005), which considered the case f'g = 
/g, a factor 2 was omitted from the Ihs of the correspond- 
ing equation (19). The subsequent algebra is nevertheless 
correct.] The assumption Ve « cr leads to a contradiction 
{ve — 0.02c[/c//g] >> (j), SO the equation has the approxi- 
mate solution 



2r7/c 2 
(J c 



-,1/3 



3/s 



\l/3 



km s 



(13) 



This solution is an attractor. At radii R large enough 
that Compton cooling becomes ineffective, the extra gas 
pressure makes the previously momentum-driven shock pat- 
tern accelerate to this value. 

At still larger radii, it may happen that the quasar sup- 
plying the driving term on the Ihs of equation (|lip switches 
off. Evidently the shock pattern will continue to propagate 
outwards for a time, because of the residual gas pressure in 
the shocked wind. Its equation of motion now becomes 







(14) 



As the independent variable t does not appear in this equa- 
tion, we let R = p, and replace R = pp', R = p^p" + pp'^ , 
where the primes denote differentiation wrt R. After a little 
algebra, the equation takes the form 



^^(pp') + iRpP + + 5cr^ 







1 r>2 / 

-^R y 



^\Ry' 



3 2 



0, 



(15) 



(16) 



where y = p^ . Now we write y = yi — 10(t^/3 to reduce the 
equation to the algebraically homogeneous form 



R^y'i + dRy'i + Gyi 







which has linearly independent solutions yi oc R 
Reversing the earlier substitutions we have 



(17) 
,R-\ 



We now choose the constants ^2,^3 to fulfil the boundary 
conditions R — 0, R = Ve at the shock position R = Rq 
where the quasar turns off. This gives finally 



if = 3 + 



10^2\ fl 2_\ _ 10^2 

3^ )\x^ 3x3 / 3 



(19) 



where x = R/Rq > 1. Figure (1) shows numerical solutions 
of the full equation of motion. With an arbitrary initial con- 
dition at small R, the shock pattern rapidly adopts the con- 
stant velocity v^. Once the quasar switches off, the velocity 
decays as predicted by the exact solution (|19[). 

Equation ((T9| gives the velocity of the shock pattern 
after the quasar switches off. This pattern stalls (i.e. R — 0) 
when 



1 



lOcr^ 



x2 3a;3 9(t;| + 10aV3) ■ 



(20) 



Since Ve » a we must have a; >> 1, so we can neglect the 
term on the rhs of pOll to get 



2 



10 V cr2 3 



9vi 
10ct2 



(21) 



where we have used Ve >> a at the last step. So finally 



Rsta 



0.95 — i?o 



: 0.95 



1 2/3 



3/gO- 



Ro 



(22) 



We can find a good approximation for the delay between 
quasar turnoff and the shock stalling by integrating eq. (|19p . 
Again neglecting the l/x'^ term this reduces to a quadrature 
of the form 

-(g/g)^^^ RdR 

[C-DK^yf^ 



(23) 



Ro 



with 



We find 



^ o / 2 10 2 \ d2 n 10 2 
C = ?,\v^ + —a jRo, D = —a 



t 



(24) 



(25) 



D 2cr2 2(7 

The shock pattern moves at the speed Ve for almost all 
the time that the quasar is on, so we can write 



Ro ~ Weia 



(26) 



where face is the timescale over which the central black hole 
accretes at the Eddington rate. Using (|22|) we can rewrite 
this as 



-Rstall 



-ta 



which of course implies 



istall — 1 



•-ace 



(27) 



(28) 



This last relation is interesting, because it shows that 
outfiows persist for quite a long time after the quasar 
switches off. Using (|13|) we find 

4taii WUcca-oTifc/fsf^^ (29) 

Hence outflows can in principle persist for an order of mag- 
nitude longer than the driving phases giving rise to them. 



5 ESCAPE 

We can use the results of the last Section to find the condi- 
tions for SMBH growth to remove gas from the host galaxy 
bulge. Attempts to explain the relation between SMBH and 
bulge mass (e.g. King, 2003; 2005; Silk & Nusser, 2010) of- 
ten invoke this kind of process. A complication so far not 
treated is that energy-driven outfiows are Rayleigh-Taylor 
unstable, and the bulge mass remaining may depend on the 
nonlinear growth of these instabilities. Nevertheless it seems 
probable that significant mass removal requires much of the 
shocked gas to escape the galaxy. 

This happens if the shock pattern reaches the galaxy's 
virial radius 



Rv 



a 

7H 



atH 

7h{z) 



(30) 



before stalling. Here H — Hoh(z), with Hq the Hubble con- 
stant, and h{z) gives the redshift dependence. Requiring 
-Rstall > Rv and using (f27)l gives 

2/3 



tacc > 1 X 10° 



??0.l/9 



2/3 



(31) 
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where 770.1 = 77/O.I. This is about twice the Salpeter 
timescale for the mass growth of the SMBH, almost indepen- 
dently of other parameters. Apparently the black hole must 
grow significantly in order to remove a significant amount 
of bulge mass. 

We may compare this accretion timescale with the time 
required for the SMBH luminosity to unbind the gas in the 
galaxy. Using eq. ([3]) with R — Rv from eq. (|30p . assuming 
that the gas binding energy is Eb ~ Ma^ and the SMBH 
energy input Ebh ~ O.OS^sI/Btvir (typical for an energy- 
driven outflow) gives 

We see that the time it takes for an Eddington-limited 
accreting SMBH to inject enough energy into the gas to 
unbind it is, in principle, shorter than the Salpeter time. 
However, crucially, this luminosity has to be communicated 
to the gas in the host galaxy. Communication via an energy- 
driven wind therefore requires an accretion timescale face 
due to the wind outflow having Ve <^ c. 

If the galaxy is inside a cluster, the outflow may reheat 

the cluster gas (King 2009). In this case, the virial radius of 

a galaxy is not well deflned, but we may consider how long it 

takes for an outflow from the central cluster galaxy to reach 

1/2 

the typical cluster cooling core radius Rcovc — 150(7 j^qqq kpc, 
where aiooo is the cluster velocity dispersion in units of 1000 
km/s. If the galaxy has a — 200 km/s and /g ~ /c, then the 
outflow cannot propagate into the intracluster medium, as 
Ve — (Jc- However, if we take the velocity dispersion of the 
surrounding material to be similar to that in a galaxy, then 
the accretion duration is 

taccc ^ 3.4 X lO^a^oV" ij) yr, (33) 

and the stalling time is tstaii,c — 3.7 x IO^ct^q yr. This 
is the timescale on which the intracluster medium is re- 
plenished by the outflow from the central galaxy, provided 
that the outflow occurs. As long as the AGN duty cycle 
of the SMBH at the centre of that galaxy is greater than 
/ > istaii,c/iH — 2.7%, the intracluster medium is continu- 
ously replenished and reheated, as the temperature of the 
gas in the snowplough phase (the outer shock) of the out- 
flow is Tout ~ 10* K, similar to the virial temperature of the 
cluster gas. 



6 VISIBILITY 

The most favourable case for viewing outflows is when each 
quasar phase is sufficiently short that the associated outfiow 
has not left the visible galaxy by the time it stalls. If for 
example we take the visible galaxy to have a size ~ 20 kpc, 
we want -Rstaii ^ 20 kpc, which by (|27|l requires 

tacc < 1.7 X W'^cr^oTifc/hr^^' yr. (34) 

Thus short growth episodes like this are most favourable for 
seeing outflows. The fraction of galaxies actually showing 
outflows then depends on the growth time of their black 
holes. The frequency of detectable outflows in principle offers 



a way of constraining the growth history of supermassive 
black holes. 



7 DISCUSSION 

This paper has discussed massive outflows in galaxy bulges, 
chiefly those driven by accretion episodes where the central 
supermassive black hole reaches the Eddington limit. We 
have shown that these outflows are momentum-driven at 
sizes R"^^ kpc, as required to explain the M — a relation, 
but become energy-driven at larger radii because the quasar 
radiation field becomes too dilute to cool the wind shock 
within the flow time. Radiative cooling is incapable of doing 
this in any regime, contrary to recent claims. 

We derive an analytic solution of the equation govern- 
ing the motion of an energy-driven shell after the central 
source has turned off. We show that the thermal energy in 
the shocked wind is able to drive further expansion for a time 
typical 10 times longer than the original driving time. Out- 
flows observed at large radii with no active central source 
probably result from an earlier short (few Myr) active phase 
of this source. 

Energy-driven outflows from longer-lasting accretion 
episodes escape the galaxy, and may well be responsible for 
removing ambient gas from the bulge, as required in some 
pictures of the black hole - bulge stellar mass relation. We 
stress however that since these outflows are Rayleigh-Taylor 
unstable, some gas make leak through the shocks and not be 
swept out. This problem is impossible to handle analytically 
and is currently numerically intractable. The inherent diffi- 
culty is that the instability sets in at very short wavelengths, 
placing great demands on spatial resolution. 
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